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ABSTRACT 



The three elastic constants of tantalum, tungsten and molybdenum 
have been measured in the temperature range between 4.2 and 300°K. The 
directly measured quantities and molybdenum and the 

computed bulk moduli of all three "elements show an anomalous temperature 
dependence. These data along with published results for other metals 
generally show a difference in temperature dependence of the bulk modulus 
between 3d and 4d-5d metals. Overlap contributions to the elastic con- 
stants are found to favor an antiferromagnetic ordering. The Debye 
characteristic temperatures calculated from elastic constants are found 
to be in substantial agreement with results obtained from calorimetry. 
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INTRODUCTION 



The elastic constants of a material determine the linear coefficients 
relating stress and strain in that material. In themselves they constitute 
an important equilibrium property. Because the equilibrium properties of 
a system are related through considerations of the system's total internal 
energy, knowledge of the elastic constants is important to a complete un- 
derstanding of thermal, electrical and magnetic behavior. 

The elastic constants of a solid are a measure of the curvature of 
the cohesive energy and their measurement yields information about the 
interatomic forces. In particular, low temperature observations are of 
importance for comparison with the results of theoretical calculations. 

Also, the low temperature elastic constants may be used to predict the 
lattice specific heat, which may then be compared with calorimetric results. 

The transition metals have always been of great interest because of 
their incomplete d shells. The Group VB and Group VIB metals, vanadium, 
chromium, niobium, molybdenum, tantalum, and tungsten, have the particular- 
ly interesting features of exhibiting a body-centered cubic structure and 
very high cohesive energies. Until recently very few measurements have 
been made on these metals, principally because of metallurgical difficulties 
in obtaining pure single crystals. 
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EXPERIMENTAL PROCEDURE 



The adiabatic elastic constants of tantalum, tungsten, and molybdenum 
have been obtained from room temperature down to liquid helium temperatures 
by measuring the transit times of elastic waves in single crystals. Un- 
rectified waveforms were observed using an ultrasonic apparatus similar 

1 2 

to that described previously.* A preliminary report on the elastic con- 
stants of tantalum pointing out anomalies in the observed bulk modulus 
3 

has already been given. 

The tantalum specimen was purchased from Materials Research Corpo- 
ration as a quarter inch diameter rod from which a 2.0643 cm (at 20°C) 
specimen was cut with accurately parallel faces. Because of ifs small 
diameter, no attempt was made to orient this crystal. Instead the orien- 
tation was taken "as is" . This orientation was determined by Laue back- 
reflection x-ray methods to deviate from about 8°, the asso- 

ciated polar angles being 0 = 84.9° and 0 = 39<.7°» Transit time measure- 
ments were made on the tantalum specimen both before and after annealing 
at 1200°C for five hours in a vacuum of at least 10 ^ mm Hg. Annealing 
did not affect the observed transit times. No change in velocity or at- 
tenuation of the elastic waves was observed in passing through the super- 
conducting transition temperature for tantalum. The transducer binder 
used with the tantalum sample was glycerin below 0°C and salol at and 
above 0°C . 

Temperatures on all samples were controlled to within 0.5°C during 
transit time measurements using a Brown Electronik recorder. A copper- 
constantan thermocouple was used for all temperature measurements above 
liquid air temperatures. For temperatures from liquid air to liquid 
helium temperatures a copper resistance thermometer was used when making 
the measurements on tantalum; a copper-gold 2.1 a/o cobalt alloy ther- 
mocouple was used for the tungsten and molybdenum measurements in the 
same low temperature range. 

The tungsten single crystal was purchased from the Linde Company 
as a 1/2" diameter rod that had been grown with its longitudinal axis 
2° off • From this a specimen was prepared that had two parallel 

faces accurately perpendicular to 2.4866 cm apart (at 0°C) . 

Nonaq stopcock grease was used as the binder for tungsten at all tem- 
peratures . 
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The molybdenum single crystal was obtained from the Wah Chang Corpo- 
ration as a rough surfaced bar averaging about one inch in diameter. Using 
standard Laue back-reflection techniques a direction was located 

in this specimen. Two parallel faces were then cut perpendicular to this 
direction, the mean distance between faces being 1.8322 cm at 0°C. The 
maximum deviation from parallelism of these surfaces in the molybdenum 
sample, as was also true for the tantalum and tungsten specimens, was 
0.0002 cm. The transducer binder used for molybdenum was Nonaq stopcock 
grease for the readings at and just below room temperature and Dow Corning 
High Vacuum Grease from below these temperatures down to helium temperatures. 

In calculating the elastic constants of these metals the thermal ex- 

4 

pension data of Nix and MacNair were used down to the limits of their va- 
lidity (the liquid air temperature region) . Extrapolations from these data 
following Grueneisen plots were used down to 4.2°K. The lattice constant 
used for tantalum was 3<,3058X at 25°C, as given in a survey by Swanson and 

Tatge. From this a density of 16.6324 grams per cubic centimeter at 0°C 

6 

was calculated. For tungsten the lattice constant given in Pearson of 
3.1650X at 25°C was used. The 0°C density of tungsten, using this value 
of lattice parameter, is 19.2710 grams per cubic centimeter. Pearson was 
also the source for the lattice constant of molybdenum, 3.1468^ at 20°C. 

The 0°K density of molybdenum calculated from this value is 10.2549 grams 
per cubic centimeter. 

Results of impurity analyses conducted on these three metals are con- 
tained in Table I. The analyses of tantalum and molybdenum were conducted 
by the firm of Kennard and Drake of Los Angeles, using test material that 
had been taken from the original specimens. The tungsten analysis, pro- 
vided by the Linde Company, is an average analysis for material produced 
at the same time as the sample used. Table I shows a low metal impurity 
content for all three samples. The interstitial gas content for tungsten 
is low, larger for tantalum and quite substantial for molybdenum. 
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TABLE lo Chemical and vacuum fusion analyses of crystal specimens 
(ppm by weight) 





TANTALUM 


TUNGSTEN 


MOLYBDENUM 



Fe 


20 


C 


10 


Fe 


50 


Cu 


3 


S 


10 


Si 


77 


°2 


222 


A1 


<100 


°2 


930 


”2 


124 


Be 


<100 


«2 


140 


«2 


9 


Cd 


<100 


«2 


10 






Mg 


<1 










Si 


<100 










°2 


25 










”2 


<20 










«2 


1 







All other impurities nil 
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EXPERIMENTAL RESULTS 



For the tantalum specimen the measured transit times were used to 
2 

calculate the quantity pV for each of the three modes being propagated 



these raw data of pV versus temperature, yielding the curves of Fig. 1. 
Values were taken from these curves at regular increments of temperature 
and a successive approximation method due to Neighbours^ was applied to 
derive the values of the elastic constants listed in Table II. 

Since the direction of elastic wave propagation in the tungsten and 
molybdenum samples was directly along > the values of calcu- 

lated from the transit time measurements are equal to certain linear com- 
binations of the three adiabatic elastic constants. Accordingly, Fig. 2 

displays the measured variation with temperature for the longitudinal mode 
2 1 

=s '2^^11'*'^12^ ^ ^44 tungsten. Fig. 3 displays the measured 
variation with temperature of the ^V^'s associated with the two trans- 
verse modes for tungsten. These equal ■2^^n“^x2^ mode polarized 

perpendicular to the z-axis and for the mode polarized parallel to 
the z-axis. Smooth curves have been drawn through the observed data 
points of FigSo 2 and 3 and a compilation is presented in Table III of 
the adiabatic elastic constants of tungsten for values taken at regular 
intervals of temperature from these smooth curves » 

Correspondingly Figs. 4 and 5 display the measured longitudinal and 
2 

transverse for molybdenum. Table IV shows the adiabatic elastic 

constants of molybdenum derived from the curves depicted in Figs. 4 and 5. 

Taking into account the various sources of error, we estimate the 

2 

absolute accuracy of a single observation of pV to be + 0.5 percent, 
although the internal consistency of the data is much better than this. 
Since the elastic wave propagation in the tungsten and molybdenum samples 
was directly along a principal direction, the above estimate of + 0.5 
percent should be valid as an estimate of error for the final reported 
values of the elastic constants. In the case of tantalum, however, the 
direction of propagation was some 8° off • Errors attendant to 

determining the direction cosines of this direction necessitate degrading 
the estimate of absolute error of the tantalum results to + 1.0 percent. 



in a direction some 8° off 




Smooth lines were drawn through 
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MEASURED ELASTIC MODULI [x 10“ dyne -cm’ 




2 

Fig. 1 Measured values of for tantalum as a function of temperature. 

The direction of elastic wave propagation for this sample was 8° 
removed from (See Text). 
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TABLE II. The elastic constants of tantalum, in units of 10 dyne-cm , 
as a function of absolute temperature. Calculated by a perturbation method 
from smooth curves drawn through values obtained from measurements made in 
a direction of propagation 8° from QllO] . The number of significant 
figures shows the internal consistency of the data and does not indicate 
the accuracy of the absolute value. 



(°K) 




«=12 


^=44 


2^^1l“^12^ 


3(Cii.2C^2) 


^12-C- 


0 


2.6632 


1.5816 


0.8736 


0.5408 


1.9421 


0.7080 


20 


2c6633 


1.5823 


0.8728 


0.5405 


1.9426 


0,7095 


40 


2.6641 


1.5844 


0.8709 


0.5398 


1.9443 


0.7135 


60 


2.6649 


1.5870 


0.8682 


0.5390 


1.9463 


0.7188 


80 


2.6652 


1.5896 


0.8647 


0.5378 


1.9481 


0.7249 


100 


2.6638 


1.5912 


0.8609 


0.5363 


1.9487 


0.7303 


120 


2.6603 


1.5909 


0.8564 


0.5347 


1.9474 


0.7345 


140 


2.6562 


1.5900 


0.8522 


0.5331 


1.9454 


0.7378 


160 


2.6518 


1.5888 


0.8480 


0.5315 


1.9431 


0.7408 


180 


2.6467 


1.5875 


0.8436 


0.5296 


1.9406 


0.7439 


200 


2.6412 


1.5860 


0.8393 


0.5276 


1.9377 


0.7467 


220 


2.6350 


1.5840 


0.8352 


0.5255 


1.9343 


0.7488 


240 


2.6286 


1.5818 


0.8309 


0.5234 


1.9307 


0.7509 


260 


2.6221 


1.5794 


0.8266 


0.5214 


1.9270 


0.7528 


280 


2.6155 


1.5768 


0.8224 


0.5194 


1.9230 


0.7544 


300 


2.6091 


1.5743 


0.8182 


0.5174 


1.9192 


0.7561 
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The elastic constant -fC^ , for tungsten as a function of temperature 
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Fig. 3 The shear constants^ C,, and 'o(C,|-C.„), for tungsten as functions of temperature. 
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TABLE III. The elastic constants of tungsten, in units of 10 dyne-cm , 
as a function of absolute temperature. The first three columns are readings 
taken from smooth curves drawn through the experimentally determined points. 
The number of significant figures shows the internal consistency of the data 
and does not indicate the absolute value. 



T 

(°K) 




^44 


1 

2(C,i-Ci2) 




"l2 


B 


^12'^44 


0 


5.3188 


1.6313 


1.6380 


5.3255 


2.0495 


3.1415 


0.4182 


20 


5.3188 


1.6313 


1.6380 


5.3255 


2.0495 


3.1415 


0.4182 


40 


5.3187 


1.6313 


1.6380 


5.3254 


2.0494 


3.1414 


0.4181 


60 


5.3187 


1.6312 


1.6379 


5.3254 


2.0496 


3.1415 


0.4184 


80 


5.3186 


1.6312 


1.6378 


5.3252 


2.0496 


3.1415 


0.4184 


100 


5.3181 


1.6310 


1.6368 


5.3239 


2.0503 


3.1415 


0.4193 


120 


5.3166 


1.6298 


1.6352 


5.3220 


2.0516 


3.1417 


0.4218 


140 


5.3138 


1.6273 


1.6315 


5.3180 


2.0550 


3.1427 


0.4277 


160 


5.3100 


1.6238 


1.6270 


5.3132 


2.0592 


3.1439 


0.4354 


180 


5.3036 


1.6207 


1.6220 


5.3049 


2.0609 


3.1422 


0.4402 


200 


5.2942 


1.6178 


1.6172 


5.2936 


2.0592 


3.1373 


0.4414 


220 


5.2846 


1.6153 


1.6124 


5.2817 


2.0569 


3.1318 


0.4413 


240 


5.2750 


1.6130 


1.6077 


5.2697 


2.0543 


3.1261 


0.4406 


260 


5.2653 


1.6108 


1.6031 


5.2576 


2.0514 


3.1201 


0.4406 


280 


5.2557 


1.6089 


1.5984 


5.2452 


2.0484 


3.1140 


0.4395 


300 


5.2462 


1.6072 


1.5937 


5.2327 


2.0453 


3.1078 


0.4381 
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Fig. 4 The elastic constant "^^44 aiolybdenuin as a function of temperature. Two separate 

measurenents of the constant, as indicated by open and closed circles, were made between liquid 
air and room temperatures . 
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Fig. 5 The shear constants, and molybdenum as functions of temperature. 
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TABLE IV. The elastic constants of molybdenum, in units of 10 dyne-cm , 
as a function of absolute temperature. The first three columns are readings 
taken from smooth curves drawn through the experimentally determined points. 
The number of significant figures shows the internal consistency of the data 
and does not indicate the absolute value. 



T 

(°K) 


"l 


=44 


l(=ll-=12> 


=11 


=12 


B 

i(=ll-2C^2> 


=12'=44 


0 


4.3650 


1.2503 


1.3855 


4.5002 


1.7292 


2.6529 


0,4789 


20 


4.3639 


1.2497 


1.3855 


4.4997 


1.7287 


2.6524 


0.4790 


40 


4.3600 


1.2488 


1.3850 


4.4962 


1.7262 


2.6495 


0.4774 


60 


4.3557 


1.2479 


1.3840 


4.4918 


1.7238 


2.6465 


0.4759 


80 


4.3548 


1.2468 


1.3820 


4.4900 


1.7260 


2.6473 


0.4792 


100 


4.3495 


1.2456 


1.3807 


4.4846 


1.7232 


2.6437 


0.4776 


120 


4.3389 


1.2446 


1.3770 


4.4713 


1.7173 


2.6353 


0.4727 


140 


4.3292 


1.2429 


1.3719 


4.4582 


1.7144 


2.6290 


0.4715 


160 


4.3212 


1.2402 


1.3672 


4.4482 


1.7138 


2.6253 


0.4736 


180 


4.3160 


1.2368 


1.3662 


4.4454 


1.7130 


2.6238 


0.4762 


200 


4^3146 


1.2333 


1.3638 


4.4451 


1.7175 


2.6267 


0.4842 


220 


4.3127 


1.2296 


1.3595 


4.4426 


1.7236 


2.6299 


0.4940 


240 


4.3061 


1.2257 


1.3540 


4.4344 


1.7264 


2.6291 


0.5007 


260 


4.2988 


1.2220 


1.3487 


4.4255 


1,7281 


2.6272 


0.5061 


280 


4.2905 


1.2188 


1.3444 


4.4161 


1.7273 


2.6236 


0.5085 


300 


4.2825 


1.2165 


1.3417 


4.4077 


1.7243 


2.6188 


0.5078 
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DISCUSSION OF RESULTS 



The elastic behavior of a cubic crystal is completely specified by 
the two shear moduli, C = and C'= reciprocal com- 
pressibility or bulk modulus, B = * C is a measure of the re- 
sistance to shear deformation across the (100) plane in the direct- 

ion and C' is the resistance to shear deformation across the (110) plane 
in the [lioj direction. B represents the resistance to deformation by a 
uniform hydrostatic pressure. 

For tantalum the computed values of C and C* are shown plotted against 
temperature as the solid lines of Figs. 6 and 7. Fig. 8 shows a corre- 
sponding solid line for B. These plots for C and C' show that these 
quantities vary with temperature in a manner similar to that observed for 
other metals. They approach 0°K with zero slope and depend linearly on 
temperature for a good portion of the temperature range. However B as a 

function of temperature (Fig. 8) shows a departure from previously ob- 

o 

served behavior by displaying a positive slope up to about 100 K, then a 
slope reversal and negative slope up to room temperature. 

In obtaining the elastic constants from a perturbation calculation it 
would be expected that, if the anisotropy is not a strong function of tem- 
perature, the perturbation corrections would be small and of the same sign 
at all temperatures. Thus if the experimental points decrease monotonically 
with increasing temperature as in Fig. 1, the elastic constants which only 
differ by small amounts from the experimental data should exhibit the same 
temperature dependence. That this is not the case for B is shown by the 
solid line in Fig. 8 which is plotted from the elastic constants listed in 
Table II. It is possible that a fortuitous choice of propagation direction 
combined with the inevitable errors in angle and transit times may lead to 
the observed behavior. 

In order to observe the effects of errors of measurement or orientation 
on the approximation method used in obtaining the elastic constants of 
tantalum, the input data were slightly changed and the elastic constants 

2 

recalculated using a GDC 1604 computer. Thus, the measured values of pV 
were individually and collectively changed by + IZ, the polar angles being 
held constant, or a measured polar angle 9 or 0 was changed by + 2°, the 
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Fig. 7 The calculated shear constant tantalum as a function of 

temperature (solid line) . The dashed lines show versus T 

when 0, or 0 are changed by the indicated amounts. 
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ig. 8 The calculated bulk modulus, tantalum as a function 

of temperature (solid line). The dashed lines show versus 

T when joV^, 0, or 0 are changed by the indicated amounts. 



values of and the other angle being held constant. The results of 

these calculations are indicated as dashed lines in Figs. 6, 7, and 8 
and show that for substantial input data changes the curves for C, C' , 
and B are merely shifted up or down their respective ordinates with no 
changes in shape. Because there is no change of shape with these per- 
turbations of the input data, it is concluded that the observed slope 
reversal effect in the bulk modulus curve for tantalum is a real effect. 

In the tungsten crystal C and C' were measured directly since pro^ 
pagation was along [noj . These results are shown in Fig. 3. The 
variation of B with temperature has been calculated and is shown as the 
upper curve in Fig. 9. 

The most striking feature of the variations of the measured quanti- 
ties C, C' and -^2) "^^44 temperature is their very flat be- 

havior from liquid air down to liquid helium temperatures. This is the 
result of there being no observable difference in the elastic wave transit 
times for the three modes from liquid air to liquid helium temperatures. 
Because of this, the bulk modulus, B, remains constant with increasing 
temperature out to approximately 100°K. Beyond this point B displays a 
positive slope out to 160°K when it takes up a negative slope and a linear 
variation with temperature. 

Fig. 3 shows tungsten is nearly isotropic in the temperature region 
measured, being exactly so at the cross-over point between the two modes 
at about 195°K. Above 140°K the shear modes vary quite linearly with 
temperature . 

As was the case with tungsten, C and C' were measured directly for 
molybdenum. These results are displayed in Fig. 5, which shows that C 
varies monotonically with temperature but C' displays a tendency to level 
off in the 90°K and 180°K regions. 

The variation of B with temperature has been calculated and is shown 
as the lower curve in Fig. 9. It is seen that the slope of the bulk modulus 
curve goes from zero to negative as the temperature rises from 0°K but there 
are two subsequent pairs of gradual slope reversals over the remainder of 
the temperature range measured, one reversal at about 80°K and one at about 
230°K. Reference to the 1% error flags in Figs. 5 and 9 readily establishes 



18 




o 

o 

ro 



O 

m 

CM 



O 

o 

CM 



o 

in 



o 

o 



o 

in 






H 



Fig. 9 The adiabatic bulk moduli "oCC, ^+20, 2) » tungsten and molybdenum as functions of temperature. 



that these slope variations of C' and B all lie within about 0.27o of a 
line that would show monotonically decreasing behavior with temperature. 
However as Fig. 4 indicates, two separate measurements of the combination 
■’■^44 versus temperature (the longitudinal acoustic mode for 
[lioj ) were made from liquid air to room temperatures, with consistent 
results. This consistency indicates that the depicted rather than the 
monotonic behavior of B versus T is warranted. 

It is seen from Fig. 1-5 that the temperature variation of the ob- 
served quantities is similar to that observed for other metals. That is 
the elastic constants approach 0°K with zero slope and, except for molybdenum, 
are remarkably linear over a wide temperature range. The percentage change 
of the constants relative to their 0°K values is, with the exception of 
tantalum, about the same as that observed for the noble metals. Near room 
temperature all the constants vary linearly with temperature and have 
significantly smaller logarithmic derivatives than the noble metals. Un- 
fortunately, the change of elastic constants with pressure have not yet 
been measured for these materials so that it is not possible to evaluate 

the intrinsic temperature dependence of the elastic constants. 

8 9 

Other investigators ’ have reported results for these metals obtained 
solely by dynamic methods. Also reported in the literature^^* are re- 
sults obtained by dynamic measurements methods combined with Bridgman's 
12 

compressibility data. 

Table V compares room temperature values from the literature with the 
results of this paper. The differences between the tantalum results re= 
ported by Bolef and by this paper tend to be close to the sum of the re- 
spective estimates of error ( 27o) . This represents fair agreement. The 
tungsten results reported by Bolef and de Klerk are in quite good agree- 
ment with those reported in this paper (within 17.) . 

We note that for tungsten the difference between Bridgman's value of 
the compressibility and the result of this work is very close to be cal- 
culated difference between the room temperature adiabatic and isothermal 
compressibilities. For tantalum and molybdenum, however, no such agree- 
ment is found, using either set of ultrasonically derived constants. 
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TABLE V. The reported room temperature values of the adiabatic elastic 2 
constants of tantalum, tungsten, and molybdenum, in units of 10^^ dyne-cm 



Element 


C 


C’ 


B 


Investigator 


Tantalum 


0.825 


0.53 


1.96 


Bolef^ 




0.818 


0.518 


1.919 


This Paper 


Tungsten 


lc51 


1.51 


2.99 


Wright^ 




1.53 


1.54 


3.08 


Bridgman^ 




1.604 


1.598 


3.083 


Bolef & de Klerk^ 




1.607 


1.594 


3.108 


This Paper 


Molybdenum 


1.09 


1.4 


2.73 


e 

Druyvesteyn 




1.068 


1.510 


2.683 


Bolef 6c de Klerk^ 




1.216 


1.342 


2.619 


This Paper 



a 


See 


Reference 


8 


b 


See 


Reference 


11 


c 


See 


Reference 


12 


d 


See 


Reference 


9 


e 


See 


Reference 


10 
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The difference in the values of the elastic shear constants reported 

for molybdenum by Bolef and de Klerk and by this paper are quite large. 

It seems clear, particularly in light of the agreement on tungsten, that 
these differences in moduli are real and that they depend on physical 
differences between the specimens rather than on differences in the mea- 

suring apparatus. That this must be the case is made more credible by 

considering that in measuring which is measured directly for [no] 

propagation, the difference in the two reported values corresponds to a 
difference in transit time of 0.5 ^s in the sample used, far outside the 
region of any reading error. 

A possible explanation of the difference in reported room temperature 
elastic constants {•^ 12 % differences for C and C' from the values of Bolef 
and de Klerk) for molybdenum is the influence of alloying of the impurities. 
Small amounts of a second phase have been found to change the single crystal 
elastic constants of copper , silver, and lithium^^. These changes 
in elastic constant have been interpreted in terras of changes in the long 
range and short range contributions to the elastic constants. The experi“ 
ments indicate that the percent change in modulus per percent change in 

composition is seldom more than four and usually near unity. Since for 

molybdenum the observed percent change in modulus per percent change in 

composition is about 100, it is difficult to see how these differences can 

be interpreted on the basis of simple alloying. 

Anelastic effects due to dislocations also contribute to change the 

moduli. A dislocation contribution (diminution) to the moduli of about one 

17 18 

percent has been observed for a single crystal of copper and of lead 
of purity in the 10 to 200 ppm range. In these experiments, the addition 

of small amounts of impurity tended to raise the observed modulus through 

19 

dislocation pinning. Chambers and Schultz have observed an internal 
friction spectrum for molybdenum in the low temperature range that shows a 
pronounced peak in the region just below 100°K (the peak shifts with 
annealing) and a high level minimum at about 200°K. They explain certain 
of the observed internal friction effects as being due to the interactions 
of gaseous interstitials and dislocations. The assay of the present moly- 
bdenum sample (Table I) shows a significantly higher gas content than the 
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sample of Bolef and de Klerk. In addition, the present sample is much 

larger, and since dislocation densities increase rapidly with the dia*= 

20 

meter of the specimen , conditions are right in the present sample for 
an anelastic mechanism such as that proposed by Chambers and Schultz to 
be quite pronounced. 

The differences in the reported elastic shear constants for moly- 
bdenum are in opposite directions, and if these differences can be ex- 
plained by means of an anelastic mechanism, then these differen^ces should 

yield some evidence as to how dislocations are distributed in the lattice. 
21 

Seeger has conducted an analysis of low temperature internal friction 

peaks wherein only a small fraction of the dislocations are considered to 

22 

lie along crystallographic directions. The work of Thompson and Holmes 
on copper indicates, however, that an appreciable fraction of the dislo- 
cations tend to lie along principal crystallographic directions. That 
the two orthogonal shear modes in the molybdenum sample are almost dia- 
metrically affected by this presumed dislocation-interstitial gas anelastic 
mechanism indicates that such a preferred alignment of dislocations may 
indeed be present. 
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BULK MODULUS 



As shown in Fig. 8 the bulk modulus of tantalum rises with increasing 
temperature above the 4,2°K value before beginning to decrease linearly. 
Observation of this anomalous behavior of B with temperature has led the 
authors to calculate the bulk moduli of the other transition elements 

whose elastic constants have been determined ultrasonically over the same 

23 

temperature range. The following data were used; Alers for vanadium; 

24 25 

Rayne and Chandrasekhar for iron; Alers, Neighbours, and Sato for 

26 27 

nickel; Rayne for palladium; Carroll for niobium; and this paper for 
tantalum, tungsten, and molybdenum. Curves presenting these results are 
shown in Figs. 10 and 11 where the bulk modulus has been normalized to its 
0°K value, denoted by B^. Fig. 10 shows B/B^ versus temperature for the 
body=centered cubic transition metals. Here the behavior of niobium, 
tantalum, tungsten, and molybdenum is significantly different from that 
displayed by vanadium and iron. Similarly in the f ace=centered cubic 

elements shown in Fig. 11 palladium has quite a different temperature de- 

26 

pendence from that displayed by nickel. Rayne has ascribed the departure 
from linearity at temperatures above 100°K of the elastic constants of 
palladium to the change with temperature of the contribution to the elastic 

constants from the unfilled d band. B/B values from the elastic con- 
28 1 ° 

stants of copper , silver and gold have been calculated and are included 
in Figo 11 to indicate the temperature dependence of the bulk moduli of 
metals that have closed-d shells. 

Taking monotonically decreasing behavior with temperature as the norm. 
Figs. 10 and 11 show that the 3d and closed-d shell metals follow the norm. 
Niobium and palladium, 4d metals, change the sign of their slope several 
times over the temperature range measured, displaying a strong positive 
slope as room temperature is approached. Tantalum (5d) shows an early 
positive slope, then tends to follow the norm above 100°K. Tungsten (5d) 
has a flat response from 0°K out to 100°K where a positive slope finally 
asserts itself. B/B^ for tungsten then parallels that of tantalum very 
closely on up to room temperature. Molybdenum (4d) generally follows the 
norm but with undulating departures. 

Figs. 10 and 11 indicate a significant difference between the behavior 
of the bulk moduli of the 3d metals, and the 4d-5d metals. Such a dif- 
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Pig. 10 The bulk nodull, normellzed to their respective 0 K values, for certain of the body-centered 
cubic transition petals as functions of temperature. Solid lines are used for 3d metals, 
dashed lines for 4d and 3d metals. 
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Fig. 11 The bulk moduli, normalized to their respective 0°K values, for certain 
of the f ace>centered cubic transition metals and for the noble metals. 



ference is in consonance with the observance of ferromagnetism in the 3d 
group but not the 4d°5d groups. 

One possible reason for this anomalous behavior of B with temperature 
that suggests itself is that these metals are involved in antiferromagnetic 
ordering in these temperature ranges, with associated anomalies in their 

thermal expansion coefficients. Evidence in support o^this from other ex° 

^ 29 

periments is scarce and difficult to assess. Bolef, de Klerk, and Brandt 
have reported anomalous thermal expansion behavior in the 180°K to 240°K 

region for vanadium. Fig. 10, however, shows that vanadium is a well*-behaved 

30 

3d metal in its B behavior. Hoare and Matthews found a maximum in the 
magnetic susceptibility of palladium between 80°K and 100°K yet Crangle and 
Smith recently looked for a specific heat A -point in this region for 

palladium, with no success. They also measured the specific heat of vana- 

32 

dium, confirming the results of Clusius, Franzosini, and Piesbergen that 
vanadium shows no anomaly between 175°K and 265°K, the region of Bolef, 
el al's anomalous thermal expansion behavior for the same metal. From the 
point of view that the body-centered cubic structure is the more natural 
vehicle for exhibiting antiferromagnetic behavior, and looking at the 
anomalous behavior of B/B^ versus T for the 4d-5d members of this group 
as possible evidence of this effect, then it is somewhat disconcerting that 
face=centered cubic palladium displays the same anomalous behavior. As the 
results of detailed specific heat, magnetic susceptibility, and thermal ex° 
pansion experiments for pure samples of these metals become available, the 
role of antiferromagnetic ordering in the transition metals may become 
clearer . 

33 

For body-centered and face^centered cubic structures de Launay has 
raised the possibility that equal to the bulk modulus of the 

electron gas of a metal. The temperature dependence of a free electron 
bulk modulus should show a steady decrease with increasing temperature as 
the lattice parameter increases. Reference to the columns in 

Tables II, III, and IV shows that tantalum, tungsten, and molybdenum like 
copper and silver exhibit just the opposite effect, an increase of ^i2“^UU 
with temperature. This reinforces de Launay' s original admonition that the 
bulk modulus as derived from his considerations is a phenomenological quantity. 
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It does not reflect on the validity of the Sommerfeld free electron model 

However correlation of C, «-C, , with the electronic bulk modulus of mono- 
12 34 

valent metals is still considered a possibility 

Excluded in the above calculations is the energy of the lowest state 

W^, Although its effect on the bulk modulus is usually taken to be small 

it is possible that this may not be the case. Further consideration of 

this contribution await a cellular method calculation of W . 

o 
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OVERLAP FORCES 



35 

Fuchs has accounted for the elastic constants of copper as being 
composed of two main parts: the electrostatic stiffnesses due to that 

part of the energy which can be regarded as the energy of an array of 
positive point charges immersed in a uniform negative sea; and the over- 
lap energy due to the overlap of the charge distributions of ions plus 

36 37 

the exchange energy. Following the lead of Fuchs and Zener , Isenberg 
has concluded that the major contribution to the elastic constants of the 
body-centered cubic transition metals is due to the overlap energy (the 
Coulomb interaction energy due to the overlap of the charge distribution 
of ions plus the exchange energy due to the same overlap). 

The relations giving the overlap contributions to the elastic con- 

g 

stants can be written as: 

N t 1 = 3/8 

N 1 = 3/^ 

N 2 = 

where N is the density of atoms, is the nearest neighbor interaction 
energy and is the next nearest neighbor interaction energy, and the 

primes represent differentiation. To a first approximation, the first 

37 

derivative of may be neglected , and the remaining derivatives are 
evaluated at the nearest of* next=nearest neighbor distance as appropriate. 

In applying the above formulas, only the part of the elastic constants due 
to overlap forces is included, i.e., the overlap contribution to the modu- 
lus is taken to be the observed modulus minus the contributions from all 
other sources. 

Using the room temperature values of the elastic constants of tungsten 
as measured by Wright^^, Isenberg determined that the second derivative of 
the next-nearest neighbor interaction, * , is positive, indicating a re- 
pulsive second-nearest neighbor interaction in tungsten. The values he 
obtained for the first and second derivatives of the nearest neighbor in- 
teraction W^’ and Wj^'' indicated a binding interaction for nearest neighbors. 
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Isenberg interpreted these results as a possible substantiation of Zener's 
36 

proposal of an antiferromagnetic arrangement for the body-centered cubic 
transition elements vanadium, chromium, niobium, molybdenum, tantalum, and 
tungsten. 



in Table VI o For all these transition elements except iron the results are 
similar to Isenberg' s. That is, the first derivative of the nearest neighbor 
interaction energy is negative and the second derivatives of both interaction 
energies are positive. 

The magnitudes of the second derivatives generally increase as one 
passes down and to the right in the periodic table. The magnitude of the 
first derivative of increases to the right in the periodic table. The 
derivatives for iron are different in sign and magnitude as might be ex° 
pected for a ferromagnetic material. 

r 

We note that if a simple exponential form = A^e i is chosen for 
both and ^2 then without neglecting W 2 * it can be shown 



where the subscript on the square bracket denotes the appropriate distance 
at which the quantity in the bracket is to be evaluated. The elastic con- 
stants used above are understood to include only that part due to overlap. 
For the transition metals considered here ^12^^ that the 

above relations predict W^ and W 2 both to be repulsive energies. This re- 
sult is in agreement with ascribing the large cohesive energy of molybdenum 

38 39 

and tungsten to the large number of electrons in the low lying d band. ’ 



Taking the published values of the 0 K elastic constants of vanadium , 
27 24 

niobium , and iron and the results of this paper these same derivatives 

11 -2 

have been calculated. The results in units of 10 dynes cm are shown 
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TABLE VI. Values of the derivatives of the nearest and next-nearest 
neighbor interaction energies attributable to overlap forces, as 
calculated from the 0°K elastic constants. Here N = the density of 
atoms, = the nearest neighbor interaction energy, = the next- 
nearest neighbor interaction energy and the primes indicate differ- 
entiation with respect to r. The subscripts 1 and 2 outside the 
brackets indicate evaluation of the enclosed quantities at the 

nearest and next-nearest neighbor distances, respectively. The value 
11 -2 

given are in units of 10 dyne-cm 



Element 


N [r W^'] 


1 


N ' 


1 . 


N [r^W^" 




Vanadium 


-1.8 




9.2 




15.4 




Niobium 


-3.2 




10.2 




20.2 




Tantalum 


-2.1 




20.5 




15.9 




Molybdenum 


-1.0 




26.2 




29.9 




Tungsten 


-0.8 




34.5 




34.5 




Iron 


+0.8 




20.2 




7.8 
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DEBYE TEMPERATURES 



Using Houston's method as outlined by Betts, Bhatia, and Wyman"^^ the 
Debye temperatures of tantalum, tungsten, and molybdenum have been calcu- 
lated from their 0°K elastic constants. Table VII presents these results 
and compares them with the results of calorimetric measurements. Although 
the values entered for the Debye temperatures of tantalum in Table VII 
differ by rather large amounts, the lattice specific heat given by 

0 = 261. 9°K is not in disagreement with the lower temperature part of the 

41 42 

specific heat data of White et al and Worley et al. 

43 

The data of Keesom and Desirant are too scattered to admit critical 
44 o 

comparison and Wolcott lists only the value 0^ = 245 K. Thus the 

situation for comparison of elastic and calorimetric Debye temperatures of 

45 

tantalum is similar to that for other metals ; i.e., where the calorimetric 
data is not too scattered, the lattice specific heat calculated from the 
elastic constants fits the low temperature part of the calorimetric data. 

Much the same remarks apply to tungsten and molybdenum. For tungsten 

46 

the data of Clusius and Franzosini is fitted at the low temperature end 

by 0 =s 384. 3°K. Clusius and Fronzosini list 0 = 380 °K and both Debye 

° 47 ° 

temperatures fit the data quite well. As in the case for tantalum Wolcott 

lists only a number and other data in the literature is too scattered for 

comparison. 

For molybdenum, the calculated value 0^ = 474»2°K gives a lattice 
specific heat in agreement with the lower temperature portion of the data 
of Clusius and Franzonini. The other published calorimetric data is not 
suitable for comparison. 
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TABLE VII. The Debye temperatures of tantalum, tungsten, and molybdenum 
as calculated from the elastic constants and as derived from low temperature 
specific heat data. 



Element 


0 

o 


Reference 


Tantalum 


261.9 


This Work 




255 


a 




246.5 


b 




245 


c 




231.0 


d 


Tungsten 


384 0 3 


This Work 




380 


e 




405 


c 


Molybdenum 


474.2 


This Work 




470 


c 




454 


e 




445 


f 



a See Reference 41 
b See Reference 43 
c See Reference 44 
d See Reference 42 
e See Reference 46 

f M. Horowitz and J.G. Daunt, Phys Rev 9^, 1099 (1953) 
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